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^ ; ABSTRACT 

Q ■ For an arbitrary quantum field in flat space with a planar boundary, an 

^ . entropy of entanglement, associated with correlations across the bound- 

ary, is present when the field is in its vacuum state. The vacuum state 
of the same quantum field appears thermal in Rindler space, with an as- 
p ■ sociated thermal entropy. We show that the density matrices describing 

1^ ■ the two situations are identical, and therefore that the two entropies are 

equal. We comment on the generality and significance of this result, and 
make use of it in analyzing the area and cutoff dependence of the entropy. 
/S ' The equivalence of the density matrices leads us to speculate that a pla- 

c^ ' nar boundary in Minkowski space has a classical entropy given by the 

Bekenstein-Hawking formula. 
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INTRODUCTION 

In a recent paper Srednicki^ considered the entropy of entanglement of a quantum field. 
Taking a free scalar field in its ground state, he calculated the density matrix which describes 
the state of the field outside an imaginary sphere when one chooses not to make measurements 
in the sphere's interior. He found numerically that the corresponding entropy scales as the 
surface area of the sphere. This entropy was originally studied by Bombelli et. al? as a 
quantum correction to the entropy of a black hole. 

't Hooft has studied a seemingly different source of quantum corrections to the entropy 
of a black hole arising from scalar fields propagating in the region just outside the horizon^. 
For very massive black holes this region may be approximated as flat Rindler space, in which 
quantum fields are thermally excited at the Hawking-Unruh temperature and carry a cor- 
responding thermal entropy. This entropy was likewise found to scale as the area of the 
horizon. 

The purpose of this comment is to point out that, for a sphere of infinite radius, the 
density matrix studied by Srednicki is equal to the thermal density matrix describing the 
fields studied by 't Hooft. (This is a special case of a more general result proven in the next 
section.) It follows that the two entropies are identical in this limit and should be thought 
of as the same quantity measured by different classes of observers. The possibility of this 
equivalence was suggested by Bombelli et. al? . Unfortunately the entropy is divergent; our 
discussion is therefore limited to its dependence on the area of the boundary and the degree 
of its divergence. Except in 1+1 dimensions, the coefficient of the divergence is non-universal. 
We will explain below that, as a result of subtleties involving regulators, Srednicki's actual 
computation^ cannot be quantitatively compared with any of 't Hooft's results."^ 

Other groups are investigating these issues and have reached similar conclusions^'^'^. 



PROOF OF EQUIVALENCE 

In this section we show that for the vacuum state of a field theory, the density matrix 
found by tracing over the fields in half of space is equal to the thermal density matrix de- 
scribing the fields in Rindler space. This result holds for a wide class of field theories; we will 
discuss its generality in the next section. 

Consider a field ^(t, x, yi, . . . , ydj^) in a {d± + 1) + 1 dimensional Minkowski space. We 
will describe the field configuration 4>{t,x,yi) as (/)ii{t,x,yi) for x > and as (/)L{t,x,yi) for 
x < 0; in subsequent expressions we suppress the dependence on the variables yi. At a fixed 
time to take the field to be in its ground state, with density matrix |0 >< 0|, and form the 
density matrix pji describing the state of the field (j)fi{to,x) for x > by tracing over all 
degrees of freedom (^L(to, x) located at a; < 0. 

< <P'r\pr\<Pr >= li^<pL < <Pl(P'r\o > < o\(Pl<Pr > (1) 

The entropy of entanglement is defined to be 

S ^ - V(/)R< (/)jl\pjll0gPR\(/)R> . (2) 

To show that this entropy is thermal, we introduce a Euclidean path integral represen- 
tation for the density matrix (1). First, we use a Euclidean functional integral to generate 
projections onto the vacuum. 

Vcpe^p - / dr dx d'^^y C 

The integral is over all fields (/)(r > 0, x) that at r = take on the value (/>^(x) for x < and 
4>'Fi{x) for a; > 0. Putting two of these functional integrals together and integrating over 4>'l{x) 
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leads to a representation of the density matrix element (1) as a single Euclidean functional 
integral on a space with a cut at the set of points {r = 0, x > 0}, subject to the boundary 
condition that (^(r = +e, a; > 0) = 4>'ii{x) just above the cut and (/)(r = — e, a; > 0) = 4>'pi{x) 
just below it, 



exp<j + / dr I dx d'^^yC 

Here e is infinitesimal and C is chosen to normalize Tr p^; = 1 . 

Now think of this path integral as generating time evolution under some Hamiltonian. In 
this geometry, instead of taking the usual Hamiltonian and time slices at fixed r, it is natural 
to use angular quantization^. Introduce the Euclidean Rindler Hamiltonian Hr which is the 
generator of rotations in the r-x plane (in Minkowski space it is the Lorentz boost generator) . 
The above path integral is then seen to be < (^^|exp {—2TrHfi) |(^^ >, so pu = e~'^'^^^ .* This 
completes the formal proof that the density matrix in half of Minkowski space is a thermal 
ensemble with respect to the Rindler Hamiltonian, with an inverse temperature (3 = 2n. ** 

DISCUSSION 

This equivalence is very satisfying: the two seemingly different density matrices are 
actually the same density matrix interpreted by different observers. An inertial observer 



* This is nothing more than the standard Euclidean demonstration that the Minkowski 

vacuum looks thermal to a Rindler observer. 
** The statement that the density matrix is thermal requires clarification. Any density 

matrix has a well defined logarithm, so p = e~^ for some H. This case is special in that H 

is the simple local operator Hr. 



who chooses not to make measurements on the field at x < obtains exactly the same 
results for any experiment as a Rindler observer who is prevented by a horizon from making 
measurements at a; < 0. In particular the two entropies are formally identical, which is a 
great technical simplification, since the entropy of entanglement defined in (2) is a difficult 
object to compute, while thermodynamics in Rindler space is relatively straightforward. 

A few comments are in order about the generality of our proof. No assumptions are made 
about the field theory, except that it must be local and must possess Poincare invariance in 
the x — t plane; the theory need not be free. The proof requires that the field be in its vacuum 
state, and does not apply if some other state is chosen. It is also specific to the geometry 
which we use; the equivalence is exact only for spaces which are divided in half by a planar 
boundary. It gives a good approximation to the density matrix which Srednicki considered^, 
for which the excluded region is a finite sphere, as long as the radius of the sphere is large. 
Finite radius corrections are especially small for the entropy of entanglement, which, as we 
will show, comes predominantly from short- wavelength modes. For example, if the theory of 
a scalar field were modified with a physical cutoff at the Planck scale, making the entropy 
finite, then, for any sphere of radius much larger than the Planck length, the entropy of 
entanglement outside the sphere would equal the thermal entropy outside a black hole of 
the same radius. Directly analyzing the entropy of entanglement in more general geometries 
requires working with the definition (2); for developments in this direction see the work of 
Callan and Wilczek^. 

To understand the physical basis for the equivalence of the two density matrices, we focus 
attention on the quadrant Q of Minkowski space with x > |t|, since Q is the part of spacetime 
causally disconnected from the t = 0, a; < half-space. Rindler and inertial observers have 



different ways of restricting their measurements to Q. A Rindler observer at a proper distance 
s from the horizon x = \t\ remains within Q by accelerating with a proper acceleration 1/s; 
consequently the observer experiences a local temperature l/2ns. An inertial observer, whose 
world line does not remain within Q, must take a different approach. Consider an experiment 
performed by an inertial observer at a distance xq from the imaginary boundary at a; = 0. 
From the moment the experiment begins, its effects begin to propagate from x = xq to x = 0; 
given sufficient time they will interact with the a; < half-space and return to a; = xq, 
allowing the x < region to affect the results of the experiment. The experiment must 
therefore be completed in a time At < 2a;o, which keeps it strictly within Q; this leads to 
energy fluctuations of order 1/2x0, comparable to the temperature experienced by Rindler 
observers at the proper distance s = xq- This highlights the connection between ordinary 
quantum fluctuations and the local temperature of Rindler space, and clarifies how the same 
density matrix can describe both situations. 

Since the density matrices are equivalent, one may compute any observable at t = 
either with the reduced density matrix (1) or with thermodynamics in Rindler space; the 
results will be identical. It follows that the half-space entropy of entanglement is formally 
equal to the Rindler space entropy. However, both entropies are divergent in most quantum 
field theories, and it seems unlikely that any observable exists which would allow them to be 
directly compared. Despite this, we will consider, for a free scalar field, the finite thermal 
entropy contained in a box in Rindler space. This quantity is adequate to study the universal 
divergence and scaling properties of the entropy of entanglement. The coefficient of the 
divergence is not universal, except in 1 + 1 dimensions, and cannot be directly compared with 
Srednicki's coefficient. 



ENTROPY CALCULATIONS 

We now present a few calculations of the Rindler entropy in a box, following the original 
work of 't Hooft"^. These results will allow us, using the equivalence proved above, to find 
the area dependence and degree of divergence of the entropy of entanglement. As we will see, 
naive scaling analysis would have given the correct answer ^'^. In more than two dimensions 
the coefficient of the divergence is non-universal. 

We begin in 1+1 dimensions. Let us compute at a fixed time the entropy contained in a 
box with walls at x = xq and x = xi;* the coordinate x measures proper distance from the 
Rindler horizon. For a well-defined counting of the number of modes within the box boundary 
conditions must be chosen at xq and xi; some of our results are sensitive to this choice. We 
work in the semiclassical (WKB) approximation; for a given (dimensionless) Rindler energy 
E the accumulated phase of the wave function is given by 

{n{E) + a)n= / —^/E^-m^x^. (3) 

Jxo ^ 

This follows from adapting equation (3.7) of 't Hooft^ for the case of large black hole mass. 
The upper limit of integration Xmax is taken to be the smaller of the large distance cutoff xi 
and the point E/m at which the integrand vanishes. The WKB quantization condition states 
that for each E such that n{E) is an integer j > 0, there is a state with energy Ej = E. The 
constant a is to be chosen according to the relevant boundary conditions. If the upper limit 
of integration is set by xi and the boundary conditions are (/){xo) = (/){xi) = 0, then a = 1; 



* Xq and xi bound the region over which a given experiment to measure the entropy is 
sensitive. They are not regulators for the field theory; in fact we perform the whole calculation 
in the continuum. 



if the upper limit of integration is set by E /m then a = 3/4 for (p{xo) = and a = 1/4 for 
4>'{xo) = 0. The WKB approximation should be quite accurate in this potential, even for low 
lying energy levels^. 

The free energy is given by a sum over energy eigenstates 

J 
the entropy in Rindler space is 5" = [/3-^ — 1 ) {l3F) with /3 set equal to 2n. For small /3mxo 
the sum over states in the partition function is dominated by highly excited states, so we may 
introduce a smoothed density of states g{E) = ^^^ ' and replace ^- — ^ j dE g{E). Note 
that we shall not see any dependence on the choice of boundary conditions, as the density of 
states is insensitive to the value of a in (3) . Integrating by parts gives 

n{E) 



(3F = -P dE 



,PE _i 



with the convention n{E) = for E < mxo- Expanding the Bose-Einstein distribution in 
powers of e~^^ and performing first the E integration, then the x integration, gives 

1 °° 1 

(3F=—y2^ [Ko{n(3mxi) - Ko{n(3mxo)] . (4) 

n=l 

In appendix A we show that these sums behave as 

7 —^Ko{nt) ~ logt + const. t <til . 

^-^ n 6 

n=l 

This gives the leading behavior of the free energy (5F in the regimes 

1 1 TT^ , Xi ,^, , 

a;o<xi<— : ---Yog V 0{mxx) 

m Tip D Xq 
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Xq <^ — <^ Xi : 
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-^log-^ + 0{1) 

D pmxo 



with corresponding entropies 

1 1 1 ^1 x/-./ 
xo <^ xi <^ — : 77 log h Olmxi) ; 

m 6 xo ._s 

111 (^^ 

Xo< — <xi: -log + C(1). 

TO D TOXo 

The first expression matches the exact 171 = result obtained in conformal field theory"^'^. In 
the second, intermediate mass case, we see that the mass takes over from xi in setting a large 
distance cutoff. 

We now consider the behavior in 1+1 dimensions for (3mxo > 1. Since (3Ej > (3mxo we 
may set 



e-OE, 



j 
In this case only the lowest lying modes in the box are occupied, and it is inappropriate to 

replace the sum over energy levels with an integral. The choice of boundary conditions is 

quite important now; for ^(xq) = 

QP ^ _g-/3mxo /g-1.84/3(mxo)i/3 _|_ ^-3.24/3(771X0)1/=^ _| \ 

while for (^'(xq) = 

QP ^ _g-/3r77Xo /g-0. 89/3(777X0)1/^ _^ g-2.59/3(r77Xo)i/^ _^ 



These energy levels are found by solving the WKB quantization condition (3) in the limit of 
large tuxq. As expected, the free energy is exponentially suppressed. If the sum over states 
were incorrectly approximated with an integral over E, the calculation would again lead to 
(4), which in the limit of large mass shows the correct e~^"^^" suppression but misses entirely 
the subleading e~^^'^^°' factors. We will see this issue is relevant for the higher dimensional 
case. 



To discuss the entropy in higher dimensions, we use a box extending from x = xq to 
X = Xi, with sides of length L in each of the d± transverse dimensions. We consider massless 
fields; finite mass effects are presented in appendix B. The free energy (3F can be expressed 
as a sum over transverse modes labeled by their transverse momentum k±. Each mode is 
equivalent to a 1+1 dimensional scalar field with an effective mass nies = \k±\, whose free 
energy (3Fi+i{meffjXQjXij (3) we have already computed. 

PF = J2 PFi+i{\k^lxo,xi,P) 

Because a long distance cutoff in the longitudinal direction is no longer needed for a finite 
result, we can set Xi = oo; in this limit, by dimensional analysis, -Fi+i is a function only of 
|/c_l|xo and (3. Approximating the sum as an integral (L/2tv)'^^ f d'^^k± and changing variables 
to K = k±xo, we have 

\2nxoJ 

This is of the form jSF — Cdj_ {/3) [^) , where Cdj_ (/3) is a dimension dependent function 
of/?. 

We now check that no important L dependence was hidden when the sum over k± was 
replaced by an integral. (By dimensional analysis this will also ensure there is no hidden xq 
dependence.) When \k\ = \k±\xo ~ xq/L, the integral should really be treated as a sum, 
whose details will depend on the boundary conditions on the walls of the transverse box. We 
can estimate the importance of this effect by computing the contribution to the free energy 
from the part of the integration region up to |k| = cxq/L, for c of order one. 

(^) " / d'^^K ^ log/3|/.| ~ log(/?cxo/L) . 



/3F = (-^] ' / d^'^K (3Fi+i{\kI(3) (6) 



Relative to the {L/xqY^ dependence of the leading term in (6), this is negligible. We therefore 
conclude that the straightforward scaling analysis is correct^'^: in any dimension d± > both 
the free energy and the entropy are proportional to the area L'^^ of the transverse box, and 
have a d±-th order power-law divergence I/xq^ as the near edge of the box approaches the 
boundary (xq -^ 0). 
The coefficient 

gets its largest contribution from modes with \k\ of order 1/13. The function /3Fi_|_i(|k| = 
meffXo,/?) increases logarithmically until \k\ ~ I//?, then decreases exponentially; additional 
powers of |k| in the integration measure push the main contribution to Cd^ (P) out somewhat 
further. As discussed above, for (3meffXo > 1 it is a poor approximation to replace the sum 
over states with an integral over a smoothed density of states. The result 't Hooft obtained^ 
for C2 (/3) may not be reliable for this reason. One can in principle determine Cdj^ (P) exactly, 
for some choice of boundary conditions, by solving the Schrodinger equation to obtain the 
spectrum of states; this probably cannot be done analytically. However, since a mode with 
large effective mass has a free energy which is quite sensitive to the boundary conditions at 
X = xo, different choices for the boundary conditions will lead to different values of Cd^(/3). 
This implies that Cd^ {(3) is not a universal quantity. 

Even if the coefficient Cd^ (P) were exactly known for a particular set of boundary condi- 
tions, it would be impossible to compare it quantitatively with Srednicki's results. Srednicki 
computes the entropy in a region xq < a; < xi with xq = a/2, where a is the spacing of the 
longitudinal lattice which he introduces as a regulator. (By "longitudinal" we refer to the 
direction perpendicular to the surface of his sphere.) In the continuum limit a — *> 0, xq — ^ 
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as well, and the entropy in his region becomes infinite. This quantity allows him to correctly 
study the scaling and the divergence structure of the entropy, but since his observable is not 
defined in the continuum, our proof cannot be used to relate his result directly to any finite 
computation in Rindler space. To see the connection with Rindler space one must define a 
measurable quantity which is finite in the continuum limit; the physical observable must be 
separated from the regulator of the field theory. 

CONCLUSIONS 

We briefiy review our main points. 

A simple proof demonstrates that the density matrix of a quantum field theory obtained 
in the vacuum state by tracing over half of space is identical to the thermal density matrix of 
the field in Rindler space. This holds for any local, renormalizable theory which is Poincare 
invariant in the x — t plane. We discuss the physical basis for this equivalence. It follows that 
the entropies studied by Srednicki^ and by 't Hooft^ are formally equal. Unfortunately, with 
few known exceptions^, these entropies are infinite, in which case it is probably impossible to 
define a physically sensible and finite observable which can be used to directly compare them. 

Using the formal equivalence of the entropies, we argue that the entropy of entanglement 
for a free scalar field is logarithmically divergent in 1 + 1 dimensions and (i_L-order divergent 
in {d± + 1) + 1 dimensions. The normalization of the entropy is universal in 1 + 1 dimensions 
and non-universal otherwise. We also show that for d± > the entropy of entanglement is 
proportional to the area of the boundary. 

We conclude by noting an important possible implication of the equivalence of the two 
density matrices. Both Srednicki and 't Hooft compute entropies associated with one-loop 
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effects of scalar fields in a background metric. In the presence of general relativity, there 
should be a tree-level contribution to both entropies. In other words, corresponding to the 
Bekenstein-Hawking entropy of the horizon of an infinitely massive black hole, there should 
be an entropy associated with an imaginary planar boundary in perfectly fiat space, with 
S = {Area)/ A. If this boundary entropy ought to be interpreted as a classical entropy of 
entanglement, then it suggests that whatever the fundamental generalization of Einstein's 
theory of gravity, it should be non-local at short distances, even at the classical level. In 
particular, Susskind and collaborators^ have suggested that in string theory this entropy is 
associated with strings which classically straddle the boundary. 

While this manuscript was in preparation, results related to those presented here ap- 
peared in refs. 4 and 6. 
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APPENDIX A: SUMS 

We consider sums of the form 






n=l 



nt] 



where Kq is a modified Bessel function. To extract the behavior for small t we calculate the 
Mellin transformation^ 



/•oo 

Jo 

= 2?-2C(2 + 
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Suppose -F(0 has poles at ^i, ^2, • • •, and write the principle part of the Laurent expansion of 
about the i^'^ pole 

Then asymptotically for small t 

fit) = E^"^^ (^1^ + ^62^(-logt) + l63,(-logt)2 + . . -V 

i 

The sum is dominated by the double pole from the gamma functions at ^ = 0, with subleading 

behavior from the simple pole in the zeta function at ^ = — 1. One finds 

2 2 

/(t)~-^logt+^(log2-7)+C'(2) + 0(t). 

For large t, f{t) is exponentially suppressed. 

APPENDIX B: FINITE MASS 

The free energy of a scalar field of mass to, 1/L <^ m <^ l/a^o? is equal to the massless 
free energy (6) plus a correction 

For (i_L < 2 the main contribution comes from the infrared, |«^| ~ toxq, where the behavior of 
the integrand is known from our 1+1 dimensional calculations. 

This gives the corrections 

^ Lm 

12— (^^ = ^^' 

-^— ^logmxo (rf± = 2). 
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Since these corrections are finite and logarithmically divergent, respectively, and come from a 
region of the integral where the integrand is insensitive to boundary conditions, we conjecture 
that their coefficients are universal. For d± > 2 the integral is dominated by \k,\ ^ mxo, 
where we may expand the integrand in powers of m. 

(J- \ d± n rl r 

= ~{mL)-'l3F\,4^-2) 
This is proportional to the free energy in dj^ — 2 dimensions, which is power-law divergent as 
xo — * and depends on boundary conditions, and so this correction is non-universal; still this 
result is rather interesting. 
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